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Using an extended Nambu-Jona-Lasinio model as a low-energy effective model of QCD, we show 
that the vacuum in a strong external magnetic field (stronger than 10^^ Tesla) experiences a sponta- 
neous phase transition to an electromagnetically superconducting state. The unexpected supercon- 
ductivity of, basically, empty space is induced by emergence of quark- ant iquark vector condensates 
with quantum numbers of electrically charged rho mesons. The superconducting phase possesses 
an anisotropic inhomogeneous structure similar to a periodic Abrikosov lattice in a type-II super- 
conductor. The superconducting vacuum is made of a new type of vortices which are topological 
defects in the charged vector condensates. The superconductivity is realized along the axis of the 
magnetic field only. We argue that this effect is absent in pure QED. 

PACS numbers: 12.38.-t, 13.40.-f, 74.90.-|-n 



Strong magnetic fields may lead to unusual effects such 
as magnetic catalysis in the (2-|-l)-dimensional Gross- 
Neveu (GN) model [Hi], in QED [3] and in QCD g]. 
The strong field supports the chiral magnetic effect in 
hot quark-giuon plasma j5] and a metalliclike conduc- 
tivity in a quarkless vacuum of lattice SU(2) Yang-Mills 
theory [6 . Recently, we suggested in Ref. [7i that an in- 
terplay between strong and electromagnetic interactions 
in a background of a sufficiently strong magnetic field 
may turn the cold vacuum into an electromagnetic super- 
conductor if the strength of the magnetic field exceeds 

Be = m^/e w lO^^Tesla, (1) 

where rup = 775.5 MeV is the mass of the p meson and e 
is the elementary electric charge. Magnetic fields of such 
strength scale may emerge in the heavy-ion collisions at 
the Large Hadron Collider in CERN ^8 and, presumably, 
in the early Universe. 

Our idea is based on a very simple argument: the parti- 
cle spectrum of QCD contains a charged vector resonance 
(a spin-triplet excitation), meson, which has a large 
magnetic dipole moment associated with an anomalous 
gyromagnetic ratio g = 2 oi the p meson. If one treats 
the p meson as a free particle, then in a background of a 
uniform magnetic field B its ground state energy [corre- 
sponding to the Lowest Landau Level (LLL)] becomes 
a decreasing function of the magnetic field strength, 
iJ^j- {B) = rn?p± — eB. The energy of p^ vanishes when 
the magnetic field reaches the value ([l]). As the field 
strength increases further, the ground state energy Ep± 
becomes purely imaginary indicating a tachyonic insta- 
bility of the ground state towards condensation of the p 
mesons. Since p^ are electrically charged, their conden- 
sation implies an electromagnetic superconductivity of 
the new ground state. Surprisingly, there is no Meissner 
effect in the g — 2 case. Moreover, the strong magnetic 
field makes the p mesons stable at QCD timescales [7]. 



The suggested vacuum superconductivity has at least 
two other analogues in particle physics: the Nielsen- 
Olesen instability of the gluonic vacuum in Yang-Mills 
theory '5], and the Ambj0rn-Olesen condensation of the 
M^-bosons induced by a strong magnetic field in the stan- 
dard electroweak model [TO] . 

In condensed matter physics, a similar phenomenon is 
known as reentrant superconductivity |llj . Usually, an 
external magnetic field suppresses the superconductivity 
via pair breaking effects, so that in a strong magnetic field 
the superconductivity is lost. However, there are super- 
conductors which may reenter the superconducting phase 
again at stronger magnetic fields. For example, in the 
uranium compound URhGe the superconductivity disap- 
pears at 2 Tesla, and then, unexpectedly, a stronger (pre- 
sumably, spin-triplet) superconductivity reappears again 
at strikingly large magnetic field of about 8 Tesla [T2]. 

There are various material-dependent proposals to de- 
scribe specific reentrant superconductors in the con- 
densed matter physics. Our suggestion in QCD [7] is 
close to the idea of Refs. [HI [13] that in a very strong 
magnetic field the Abrikosov flux lattice of a type-II su- 
perconductor may exhibit a "reentrant" quantum regime, 
characterized by the LLL dominance, the absence of the 
Meissner effect, a spin-triplet pairing, and a supercon- 
ducting current flow along the magnetic field axis only. 

In Ref. f7| we suggested the existence of the new super- 
conducting phase using effective bosonic electrodynamics 
of the p mesons of Ref. [M]. However, this model treats 
the p mesons as pointlike particles, and thus it may be 
inapplicable at strong magnetic fields ([!]) when the mag- 
netic length becomes of the order of the size of the p 
meson. Here we use the much more general fermionic 
Nambu-Jona-Lasinio (NJL) model [15] as a low-energy 
effective theory of QCD in order to show the existence 
of the electromagnetic superconductivity induced by the 
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strong magnetic fields in the vacuum. 

We consider an extended two-flavor {Nf = 2) Nambu- 
Jona-Lasinio model with three colors {Nc — 3) [16]: 

^P)^i){i$ + QA- M°)V' + + C'^^ , (2) 

where the light quarks are represented by the doublet 
tp — {u,d)^ and M'^ = diag(m° , to°) is the correspond- 
ing bare mass matrix. The uniform magnetic field back- 
ground B = (0, 0, B) is encoded in the Abelian gauge field 
A^" = iA°,A) = (0,-5x2/2,5x1/2,0), and the electric 
charges of the quarks, qu = +2e/3 and qd = ^ e/3, are 
combined into the matrix Q = diag(g„, qd)- The hat over 
a symbol indicates a 2 x 2 matrix in the flavor space. 

The last two terms in Eq. ^ represent the scalar and 
vector four-quark interactions, respectively: 



.(4) 



where G^g ^ and G 
and T — {t^, ' 

We follow the standard approach [TFI and introduce 
the following bosonic fields corresponding to the quark- 
antiquark bilinears: one scalar field cr ~ the triplet 
of three pseudoscalar fields tt ~ -ipj^filj [made of the elec- 
trically neutral, 7r° = tt'^, and electrically charged, tt* — 
(tt^ =F in^)/^/2, pions], four vector fields ~ 'ip^f^r'^'ip, 
and four pseudovector (axial) fields A^^ ~ {jj^^^^r'^ij), 
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(Vi7^rV) 



(3) 



(i^7^75rV)' ,(4) 



-y' are corresponding bare couplings, 
^) are the Pauli matrices 



V,., 



=0 ^ 



(5) 
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The vector-meson matrix ([s]) is composed of the singlet 
(in the fiavor space) vector (in the coordinate space) uj- 
meson field w^, while p° = pi and yO± = (p^ Tipl)/V^ 
represent, respectively, electrically neutral and charged 
components of the p-meson triplet. The light axial 
mesons are encoded in the matrix ([6]): the fields and 
(a°,aj) represent, respectively, the singlet axial /i me- 
son and the ai triplet of the axial mesons, respectively. 

We rewrite the four-quark interactions ([3| and Q via 
Gaussian integrals over the bosonic fields cr, tt, V^t, A^, 
and integrate over the quarks in the partition function: 



i J df^x C _ 



DaDnDVDAe'^^''-^'^-^^ 



where the effective bosonic action is 
S[a,TT,V,A] = S^[a,TT,V,A] 



(7) 
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^ -iNcTiLniiV) , (8) 
iV = i$ + QA-M° + 7^4 - (ct + «7^7?t) . (9) 



Next, we calculate the effective action ([t]) in the strong 
magnetic field background in the mean field approach. 
We use simplified notations for the expectation values of 
the fields, (a) = a etc. In the absence of the external 
magnetic field the expectation values of the fields tt, V, 
and A are zero [TB] , while the expectation value of a plays 
a role of the constituent quark mass, m.q — a ^ 300 MeV. 

In order to simplify our calculations, we notice that 
the presence of the external magnetic field breaks the 
fiavor symmetry down to its diagonal subgroup, so that 
the diagonal chiral rotations f2 = e*"^"^ can still be 
used to eliminate the neutral pion condensate tt". We 
also neglect the mass matrix because our calculations 
show that the conducting properties of the vacuum are 
almost independent of the hare quark masses d^^- 

The operator ([o]) can be represented as the sum iD = 
Wq + W of the tree-level operator iVo = i0 + Q Jl — a 
and the contribution W from the "exotic" condensates. 



W = f + 7^4 - i75(7rVi + ttV^) . 



(10) 



At low magnetic fields W = (W) = 0. Let us assume 
that, at certain strong magnetic field B ~ B^'^^, the ex- 



pectation value of the condensate ( 10 1 is nonzero. Let us 
advance slightly into the new phase taking B > B^'^^, 
so that the magnitude of the suspected condensate is 
still small, < \W\ < cr. Then the effective action ^ 
can be expanded in powers of the W field, and the fact 
of the emergence of the new condensate should be seen 
as a tachyonic instability of the corresponding potential 
at W" = (alternatively, we could also solve Schwinger- 
Dyson gap equations and get the same result). 

The tree-level propagator 5"^°) = V^^ of the fermion 
doublet in the strong magnetic field has the following 
form: 5(°)(x,y) = diag[S^°\x , y) , sf'> {x , y)] , where Sf 
is the propagator of the quark species. 

The p-meson condensation and, consequently, the in- 
duced superconductivity are the LLL phenomena [7]. 
Thus, it is natural to restrict ourselves to the LLL ap- 
proximation which usually gives a dominant contribution 
to nonperturbative low-energy quantities in the limit of 
the strong magnetic field (in particular, because of the 
large energy gap SE ^ {eBy^^ between the LLL and 
higher levels) [3l lU [17]. The full solution beyond the 
LLL approximation will be presented elsewhere. 

In the LLL regime the propagator S^^'^ factorizes into 
the _B-transverse and _B-longitudinal parts which depend, 
separately, on the i3-transverse 
longitudinal, x"=(x",x'^), coordinates [Sj: 

4°)'^^^x,2/) -P/(x^,y^)5j(xll -yll) (11) 
(below we omit the superscripts "(0)" and "LLL"). Here 
QfB\ 



(x\ x^), and B- 



27r 



-62 



(12) 



is the transverse projector onto the LLL states and qf is 
the electric charge of the f'^^ quark. 



3 



The longitudinal part of the fermion propagator (111, 
^sgn{qfB) basically, a fermion propagator in the 
1 + 1 dimensions (we always take eB > for definiteness) , 

„ ]i-z/y^2 



7II fell 



p 



f 



p 



f 



,(13) 



and the matrix Pj (we use f — ±1 for, respectively, 
f = u,d) is the spin projector operator onto the fermion 
states with the spin polarized along (for u quarks) or op- 
posite (for d quarks) to the magnetic field. The operator 

pj' projects the original four 3+1 fermionic states onto 
two (l+l)-dimensional fermionic states, so that fermions 
can move only along the axis of the magnetic field. The 
projector ( 12 ) satisfies the relation P^oPj- = P^, where 
"o" is the convolution operator in the i?-transverse space. 

For a coordinate-independent condensate cr, the zero- 
order (in powers of W) contribution to the effective 
action M gives us the potential V{a) = V^^\a) + 
a"^ /{2G g ) related to the action as 5 = — / d^x V with 



\eB\N, 



0- In ^ - 






L ^2 







where 1/e = l/e - je + log47r, je ~ 0.57722 is Euler's 
constant, and /.t is a renormalization mass scale. In order 
to regularize the divergent contributions of the (1 + 1)- 
dimensional fermions we implemented the dimensional 
regularization in d = 2 — 2e dimensions. The renormal- 
ization of the NJL coupling constant in the MS scheme, 
l/Gs = -Nc\eB\/{4Tr'^e), resembles the renormal- 

ization of the 1+1 dimensional GN model [T] with the 
identification Gs = 27rG'GN/(-^c|ei3|) first noticed in [3]. 
The minimum a = amin of the renormalized potential. 
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In 



- 1 



\eB\N, 
87r2 \ fi- 

provides us with the i?-dependent quark mass 

m,(B) - a^i„(B) = Aicxp{-27rV(GsA^c|eS|)}.(14) 

In the LLL approximation to the NJL model the scale fi 
is not fixed as it is related to the B-longitudinal 1 + 1 
motion of the quarks. Beyond the LLL approach the 
scale may perhaps be set as oc |ei?| following Ref. [1]. 

The effective bosonic model T suggests that the pos- 
sible superconducting ground state should exhibit an in- 
homogeneous behavior in the i?-transverse plane. Thus, 
we assume that the exotic condensates may be x-^- 
dependent, W = W{x-^), and calculate the correspond- 
ing quadratic contribution to the effective action (|8|, 

(15) 



c(2) _ 



W W . 



(despite the expected ^-''-dependence of the condensates 
W, we still call the functional Vl as the "potential"). 



We find that the potential ( 15 1 involves only the B 



transverse components of the vector and axial mesons, 
mJeBl 



97r2 



(J 



In 



o P„ o , 



i-ln^-2)(r °Pe°0 



(16) 



where </> = (p+ + ip^)/2 and ^ = {a^ + ia^)/2. The 
_B-transverse projector for the unit charged particle, 

P^{x^,V^) = {9n/\eB\)P.^{x^,y^)Pi{y^,x^) is given 



by Eq. (12 1 with the replacement qf — > e. 

The unstable tachyonic mode of the potential ([7|, 



( 16 1 turns out to be an inhomogeneous eigenstate of the 



charge-1 projection operator Pg, 

(Peoc/.)(a;^)==0(x^), (17) 
The solution is a general Abrikosov-like configuration [18] 

(j)^cj)oK{z/LB) , Lb = ^2^IW\. (18) 



-|-CX3 



-7rn -|-27rn2: 



(19) 



where 0o and c„ are arbitrary complex parameters and 
x^ + ix^ (and similarly for the axial vector field i^). 



z 



The solution (18) represents a (periodic) flux-tube 
structure similar to the Abrikosov lattice which is realized 
in a mixed state of a type-II superconductor subjected to 
a near-critical external magnetic field ^ISj . Generally, the 
coefficients c„ can be fine-tuned by a complicated min- 
imization procedure if the full potential is known |18) . 
Here we follow Ref. [71 [10] and set c„ = 1 so that the so- 
lution ( 18 1 represents a square lattice with the quantized 
area 27r/|eP| = L\ given by the magnetic length Lb- 



The quadratic potential, evaluated at the solution ( 18 1, 

27VJeP 



1 



(I0ol^ + ieo|^) 



97r2 



is unstable towards a spontaneous creation of the B- 
transverse condensates with the tachyonic mode p\ — 
) if the strength of the magnetic field exceeds 



B 



NJL 



9^2 



1 

Gb 



1 



(20) 



2eNcGB Gb Gy 9Gs 

with 1/Gy = l/G|^^-A^c|eP|/(97r2e') and 1/e' = 1/e-l. 

Due to the fact that the phenomenological values of the 
parameters Gsy vary in a broad region }19j . and due to 
subtleties of the renormalization of the effective (1 + 1)- 
dimensional theory embedded in 3+1 dimensions, we can 
only give an approximate estimation of the critical field: 
ePc - 1 GeV^ or B^ - lO^^ ^ggj^^ 

The quartic correction to the potential in Eq. ([8|, 



(4) 



27r2m2 

allows us to find the condensate at B ^ B^'^^: 

,1/2 



fo 



'B) = e''«C^m,{B){l-B^''^/B 



(21) 



(22) 
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where is a constant phase, Co ~ 1.2, w 0.51 and 
the quark mass to, is given in Eq. (14). At i? < Bf'^^ the 
condensate (22) is zero. The phase transition at B ~ Be 



is of the second order with the critical exponent 1/2. 
Thus, the magnetic field induces the quark condensate 



(U7irf) = -^{ui2d) = Po{B)k{^^^^) ^ p{x^), (23) 



where po{B) — (j)o{B)/Gv- Using known (see, e.g., Ref. 
[I8] ) general properties of the function K{z), Eq. (19), 
we conclude that the ground state should be given by a 
periodic (in general) lattice of a new type of topological 
vortices which are parallel to the magnetic field. The 



phase of the condensate ( 23 ) winds around the center of 
each vortex where the absolute value of p{x-^) vanishes. 

The condensate (23) locks the local U{l)c.m. trans- 
formations with the global O(2)].ot rotations of the co- 
ordinate space about the magnetic field axis [71 ED] : 
t^(l)o.m. X 0(2)i.ot Giat, where Giat is a discrete sym- 
metry group of rotations of the p-vortex lattice. 

The new vacuum state is superconducting. One can 
show that there is no i?-transverse current, J^—J^=0, 
so that the electric current flows along the magnetic 
field axis only. In a very weak (test) electric field 
E = {0,0, E^) with Ez <C B, the induced electric cur- 



rent in the new vacuum state (23) in a linear-response 



approximation is (we use the retarded Green functions): 
"^''(^) = E, Qf{^fr^f)^-tr[rQSix,x)], (24) 



We average the current (24) over the _B-transverse 



plane and, in the leading order in powers of p, we get: 



5Q(J) ^ dj{x\\) 



2C, 



dz 



dt 



(27r)= 



" eh 



B - B. 



NJL\ 



E, 



(25) 



where Q is the plane- averaged electric charge density J'^, 
J is the plane-averaged current J^, and « 1 [21]. At 
B < Be the right hand side of Eq. (25) is zero. Apart 



from prefactors, the transport laws in the NJL model ( 25 ) 



and in the p-meson electrodynamics [7] are identical. 
The linear- response law ( 25 ) can be rewritten in a 

{F,F)Fi"', via the 
E^). 



Lorentz-covariant form, d^^'J'^^ — 7 
invariants {F,F) = 4(5, i?) and {F,F) ^ 2(5^ 



Here = e^.uapF'^P /2 and 7 is a function of {F, F) [20] . 



Equation ( 25 ) is a London equation for an anisotropic 



superconductivity. Thus, we have just shown that the 
strong magnetic field induces the new electromagneti- 
cally superconducting phase of the vacuum if _B > B^. 
An empty space becomes an anisotropic superconductor. 

The superconductivity of the vacuum is a new effect 
which is realized at the QCD-QED interface. This mech- 
anism should not work in the pure QED since electrically 
charged spin-1 bound states are absent there. 

On general grounds one can expect that increase in 
temperature T (which, in general, should be of a hadronic 



scale) should lead to an evaporation of the p condensate 
with a loss of the superconductivity. The suggested low- 
T part of the B — T phase diagram is shown in Fig. [ij 



T 



Hadronic phase 



tran 



sition 



electromagnetic 
superconductivity 



B 



FIG. 1. Low-temperature part of the QCD phase diagram. 
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